Artificial kagome arrays of nanomagnets: a frozen dipolar spin ice 
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Magnetic frustration effects in artificial kagome arrays of nanomagnets are investigated using 
X-PEEM microscopy and monte carlo simulations. Spin configurations of demagnetized networks 
reveal unambiguous signatures of long range, dipolar interaction between the nanomagnets. As soon 
as the system enters the spin ice manifold, the kagome dipolar spin ice model captures the observed 
physics, while the short range kagome spin ice model fails. 

PACS numbers: 75.10.Hk, 75.50.Lk, 75.70.Cn, 75.60. Jk 



Complex architectures of nanostructures are nowa- 
days routinely elaborated using bottom- up and/or micro- 
fabrication processes. This technological capability al- 
lows scientists to engineer materials with properties that 
do not exist in nature, such as negative refractive index 
[l] or tunable photonic bandgap [2 , but also to elabo- 
rate model systems to explore fundamental issues in con- 
densed matter physics. Frustrated spin models are one 
example for which theoretical predictions can be tested 
experimentally [3]. In these models, spins are placed 
on an infinite network of a given geometry, and inter- 
act through a short or long range coupling of ferromag- 
netic or antiferromagnetic nature. Under certain condi- 
tions, magnetic frustration effects occur, giving rise to 
an infinite degeneracy of the ground state, i.e. a finite 
entropy at zero temperature. Several experimental stud- 
ies on two-dimensional, artificial arrays of ferromagnetic 
nanomagnets have been reported recently [5-5 , bringing 
together the communities of magnetic frustration, sta- 
tistical physics and nanomagnetism. Properties of these 
arrays seem to be well described by short range spin ice 
(SI) models: they show disorder, satisfy ice-like rules, 
and correlation coefficients fit with predictions [6 . 

In the kagome SI model, spins interact through first 
neighbor exchange coupling only. Experimentally how- 
ever, nanomagnets are coupled via the magnetostatic 
interaction, which extends beyond nearest neighbors. 
While the main interest for frustrated compounds arises 
from the massive degeneracy of their ground state, this 
degeneracy is fully lifted when long range, dipolar inter- 
actions are included in the model [H |10] . Understanding 
whether or not these long range interactions influence the 
local spin configuration in artificial arrays of nanomag- 
nets is thus essential, especially because these networks 
are often considered as a playground to study magnetic 
frustration effects on a mesoscopic scale. Combining ex- 
perimental and theoretical investigations, we show that 
the kagome dipolar spin ice (DSI) model captures the 
physics we observe, where the kagome SI model fails. 




FIG. 1: (Color online) (a) 15 /xm X-PEEM magnetic image 
of a demagnetized network. Direction of magnetization is de- 
duced from the black and white contrast, as illustrated for two 
vertices that satisfy the ice rule, i.e. the local lin/2out (mag- 
netic charge Q=— 1) or 2in/lout (magnetic charge Q=+l) 
spin configuration, (b) Sketch of the array indicating the rel- 
ative indices used for the 7 nearest neighbors. The partly 
superimposed dots highlight the nodes of the kagome lattice. 



Honeycomb arrays of nanomagnets have been fab- 
ricated from a Si//Ta(5nm)/Co(10nm)/Ru(2nm) thin 
film. An aluminium mask, defined by e-beam lithogra- 
phy and lift-off, has been used to protect the structures 
during ion beam etching. The Al mask has been subse- 
quently removed by chemical etching. The nanomagnets 
have a typical size of 470 x 70 x 10 nm^, ensuring a high 
uniaxial shape anisotropy, so that they can be consid- 
ered, locally, as Ising pseudo spins. The distance between 
nearest nodes in the network is set to 500 nm (Fig. la). 
The Co nanomagnets are thus physically disconnected, 
and are coupled through the magnetostatic interaction. 
The networks of nodes have kagome topology (Fig. lb) 
and are composed of 342 Co nanomagnets. The magnetic 
configuration of each element is resolved using the X-ray 
PhotoEmission Electron Microscope (X-PEEM, SOLEIL 
instrument) of the Nanospectroscopy beamline at the 
ELETTRA synchrotron facility (Trieste, Italy). Mag- 
netic images are acquired at the Co L3 edge at rema- 
nence, after demagnetizing the networks ex situ. The 
importance of the demagnetizing procedure to approach 
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TABLE I: Theoretical (SI model) and experimental spin-spin 
correlators, number of ij interactions in our networks, and 
coupling coefficients in arbitrary unit deduced from the point- 
dipole approximation and from micromagnetic calculations. 
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the ground state spin configurations has been highlighted 
in previous works [TT',n!2]. To demagnetize our networks, 
the sample is rotated with a period of about 20 ms in a 
damped, in-plane sinusoidal magnetic field (1 s period, 
0.11 Oe amplitude variation between two periods). From 
the X-PEEM images, the magnetic state of each individ- 
ual nanomagnet is measured (Fig. la). Data presented 
here have been obtained from 9 nominally identical net- 
works, which have been demagnetized 6 times, out of 
which 34 magnetic configurations have been determined. 

The kagome ice rule, i.e. the preferential selection of 
lin/2out or 2in/lout low-energy spin configurations (see 
Fig. la), is globally well obeyed: over all magnetic con- 
figurations, less than 0.6% of the vertices have a 3in or 
3out high-energy spin configuration. These "forbidden" 
vertices are likely due to structural defects in the net- 
works and often appear at the same locations. Residual 
magnetization after demagnetization is low, of the order 
of a few percents, consistent with refs. ITTIand 12, From 
the analysis of our magnetic images, we deduce the spin- 
spin correlation coefficients Cij = (Si.Sj) between nor- 
malized Si and Sj spins (indices are relative, identical to 
those used in previous literature [13] , and defined in 
Fig. lb). The experimental values, averaged over the 34 
demagnetized networks, are reported in Table [l| for the 
7 nearest neighbors. These values are compared to the 
coefficients expected from the SI model, described by the 
Hamiltonian Hsi = —Ja(3 ^{ij) ^i-^j that couples near- 
est neighbors only, through the Jaf3 coefficient [131 IE] • 
Although differences are visible, as previously observed 
experimental results and theoretical predictions from 
the SI model agree relatively well at first sight. 



FIG. 2: (Color online) (a) Experimental histogram deduced 
from the 34 demagnetized networks. Color code is used to sort 
the configurations in 3 subsets. Curves show the normalized 
^long distributions from the SI model (dotted line) and after 
replacing the expected Cau and Car] coefficients by their mea- 
sured values (continuous line), (b) Experimental correlation 
coefficients averaged within each energy subset. 
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"j, and D = (/io/47r)/i^ the strengh of 
the dipolar coupling. The corresponding Jij values are 
reported in Table |T] and compared to the values obtained 
using micromagnetic simulations . Good agreement is 
found between the two approaches, except for first neigh- 
bor interactions, for which the minimal distance between 
the nanomagnets is small compared to their extension. 

The system energy can be split in two terms: E = 
^short + ^long whcrc £^short IS the energy related to 
the interaction between first neighbors only (£^short = 
—Naf3Ja(3Ca^/N), and ^long the energy related to the 
interactions beyond the nearest neighbors. We will show 
that, once the system enters the spin ice manifold, Ca^s 
is configuration-independent, i.e ^short is constant. In 
the following, we thus focus our discussion on the £^iong 
term. The presence of a few "forbidden" vertices in our 
networks is then not critical as ^long excludes energy 
contributions from nearest neighbors. The experimen- 
tal ^long distribution is shown in Fig. 2a, and compared 
to the distribution expected from the SI model (dotted 
line). Interestingly, this distribution is shifted towards 
negative energy. To understand where energy is gained 
in our networks, we sort the 34 measured configurations 
in 3 subsets corresponding to the lowest (red) , intermedi- 
ate (green) and highest (blue) energy configurations (Fig. 
Considering the nanomagnets as Ising spins [T5], we 2a), and we average the Cij coefficients within these 3 



can derive from the Cij values the network energy E = 
— (1/A^) "^Zi^j NijJijCij. In this expression, is the to- 
tal number of spins, Nij and Jij being respectively the 
number of interactions and coupling coefficients of type 
ij (Fig. lb). Within the point-dipole approximation, 
the coupling coefficients Jij can be calculated: for two 



subsets (Fig. 2b). Two main features are then observed: 
i) independently of the configuration energy, the C^v and 
Car] coefficients have higher values than expected (also 
see Table 1), and ii) the values of the Cq,^ and C^s co- 
efficients strongly vary with energy, contrary to what is 
observed for the 5 other measured coefficients. 
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As depicted in Fig. lb, the Cocv and Cq,^ coefficients 
reflect correlations along a line in the network. We might 
think that these higher correlation coefficients originate 
from the demagnetization procedure that may favor fer- 
romagnetically coupled nanomagnets along these lines. 
However, although the system has threefold symmetry, 
the three main directions of lines are not equivalent, and 
we observe anisotropy in the Cav and Car] values (one 
direction is different from the two others). We thus at- 
tribute the high C^v and Car] values to variations of 
coupling coefficients between flrst neighbors, presumably 
due to anisotropy in the fabrication process. Replacing 
the theoretical correlation coefficients along a line by the 
experimental ones, keeping all the other coefficients un- 
changed, leads to a larger shift of the distribution towards 
negative energies (continuous line in Fig. 2a). This new 
energy shift reproduces well our data. Thus, the total 
energy shift we observe experimentally has two main ori- 
gins: an intrinsic shift associated with the ice model itself 
(^long is lower in the SI manifold than in the paramag- 
netic state in which all spin conflgurations are equiprob- 
able [17 ), and an extrinsic shift due to the lithography 
process. Both phenomena are related to flrst neighbor 
interactions, and do not involve long range coupling. 

In contrast with the other correlation coefficients, the 
values of C^^ and C^s strongly vary with £^iong (see Figs. 
2b). For several conflgurations, C^-f and C^s reach ab- 
solute values that can be surprisingly large compared to 
the values expected from the kagome SI model. Mag- 
netostatic coupling between the (a, 7) and (a, 5) spins 
favors antiparallel alignment of the corresponding nano- 
magnets, and thus negative values of their spin-spin cor- 
relators. Intuitively, this could explain why the measured 
Ca-f and Ccx5 values are not in quantitative agreement 
with the values expected from the kagome SI model. 

To test whether this interpretation better describes 
our flndings, monte carlo simulations of the kagome DSI 
model have been performed, using the following Hamil- 
tonian: i^DSi 



T^i^^Jij Si'Sj^ where 
Ji is an additional isotropic nearest neighbor coupling. 
This Ji term can be tuned so that the effective near- 
est neighbor interaction Ji + Ja/s = ^Ja(3i to account 
for the spatial extension of the nanomagnets (see Table 
1). The simulations were done for L x L x 3 sites lat- 
tices (L ranging from 11 to 144) with periodic bound- 
ary conditions, using a single spin flip algorithm and 
a simulated annealing procedure from T/J^^ = 100 to 
T/Ja(3 = 0.1. In these simulations, 10^ modifled monte 
carlo steps (mmcs) are used for thermalization, where 
1 mmcs corresponds to a set of local updates sufficiently 
long to achieve stochastic decorrelation, itself determined 
using the on the fly calculated spin-spin autocorrelation 
time. Measurements follow the thermalization process 
and are computed with 10^ mmcs. Results of the simula- 
tion are shown in Fig. 3 where spin-spin correlators are 
plotted as a function of temperature. In the high temper- 
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FIG. 3: (Color online) Correlation coefficients Cd derived 
from the DSI model as a function of temperature T. Greek 
index i — f3^^^ . . . identifies each correlator. Inset : tempera- 
ture dependance of correlators within the SI model. 



ature regime, the system is in a paramagnetic state. As 
temperature is lowered, the system enters the short range 
kagome spin ice manifold {Ca(3 = 1/6), and the 7 mea- 
sured correlation coefficients nearly reach the values ex- 
pected from the SI model. We emphasize that in absence 
of long range interaction, these values would then remain 
temperature-independent (see inset of Fig. 3). The situa- 
tion is different when dipolar interactions are included in 
the calculation. At very low temperatures, a transition 
to an ordered phase takes place [9l [10]: all correlation 
coefficients vary abruptly (not shown here), and strongly 
deviate from the SI values. In the intermediate regime, 
i.e when 2 < T/ J^p < 5 (see Fig. |3|, mean values of the 
correlation functions do not allow easy discrimination be- 
tween the SI and DSI models. However, two correlation 
coefficients slowly decrease when lowering the tempera- 
ture (see the two lowest curves). These two coefficients, 
Ca-f and Ca5, are precisely those for which we observe a 
change experimentally, when £^iong decreases. This strik- 
ingly contrasts with predictions from the SI model and 
fits well our experimental data. 

An alternative description of these spin models is use- 
ful. Identifying each spin Si to a magnetic dipole allows 
the mapping of any spin conflguration on the kagome 
lattice onto a classical magnetic charge conflguration on 
the honeycomb lattice (Fig. la). The high tempera- 
ture paramagnetic phase of the SI and DSI models then 
translates into a paramagnetic charge state. When tem- 
perature is lowered, both models enter the spin ice man- 
ifold: the kagome ice rules are satisfled, and charge val- 
ues are constrained to ±1 (instead of ±1 and ±3 at high 
temperature), therefore corresponding to an emergent Z2 
fleld. In the SI model, this Ising fleld remains weakly cor- 
related and temperature independent. The DSI model 
however behaves differently: at lower temperatures, it 
leaves a flrst spin ice manifold (SI-I) and enters a more 
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constrained phase (SI-II), in which spins stih fluctuate 
but where charges crystahize in an antiferromagnetic pat- 
tern [To]. During the phase transition of the hidden Z2 
fleld, the stochastic dynamics is strongly modifled, with 
exponentiahy suppressed single spin flips, leaving only 
loop-like moves to operate efliciently in the SI-II man- 
ifold. We have thus investigated the temperature de- 
pendance of the short range magnetic charge correlation 
{Qi'Qi+i)-, and compare monte carlo predictions from 
the SI and DSI models (Fig. 4a) to our measurements 
(Fig. 4b). In the highest temperature regime, where 
all charge states are accessible (Q=±l and Q=±3), the 
unconstrained paramagnetic charge phase is associated 
with an asymptotic correlator (Qi.Qi+i) =-1 (Note that 
a Q=±3 charge state corresponds to a monopole defect, 
see for example refs. [TSl and 119]) . When the temperature 
is lowered, the SI model enters the spin ice manifold, with 
constrained Q = ±1 charges and a temperature indepen- 
dant short range charge correlation (Qi.Qi+i)=-l/9. In 
the DSI model however, the charge correlator flrst in- 
creases to reach a maximum value around -0.15, and then 
steadily decreases with temperature. Note that there ex- 
ists a temperature range within the spin ice manifold, in 
which the mean values of the correlator, as well as their 
associated standard deviations, are mutually exclusive 
for the two models (Fig. 4a). Analysis of our X-PEEM 
images reveals strong deviations from the -1/9 value of 
the SI model (Fig. 4b), clearly outside standard devia- 
tions of the correlator distributions, and agrees well with 
predictions from the DSI model. We emphasize that in 
all 34 networks, the spin conflgurations are charge neu- 
tral, which rules out a global extrinsic breaking of the 
Z2 symmetry. Therefore, we attribute the deviations of 
the short range charge correlations to the proximity of 
the SI-II regime, in which the Z2 and translation sym- 
metries are broken. Unambiguously, our results reveal 
that dipolar interactions are at work in our arrays. 

In summary, numerical computation and magnetic 
imaging show that the kagome SI model fails to interpret 
our experiments, while the kagome DSI model captures 
the physics we observe. The fact that artiflcial arrays of 
nanomagnets may look like short range kagome spin ice 
systems at flrst sight is fortuitous, and related to the crit- 
ical slowing down at low temperature of the single spin 
flip dynamics. Therefore, while artiflcial kagome arrays 
of nanomagnets should not exhibit flnite entropy at zero 
temperature [9^ , they are most likely described by a 
frozen dipolar spin ice regime with an emergent symme- 
try breaking of the hidden Z2 fleld associated with the 
classical magnetic charges. 
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FIG. 4: (a) Theoretical short range magnetic charge cor- 
relators {Qi.Qi-\-i) (bold line) and their standard deviations 
(hatched regions) calculated from the SI and DSI models as a 
function of temperature, (b) Experimental {Qi.Qi-\-i) values 
measured for the 34 demagnetized networks as a function of 
^long- (c) Histogram of the experimental charge correlators 
compared to the nearly gaussian distribution computed in the 
SI model, within the spin ice manifold. 
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